Relaxation and dephasing in open quantum systems time-dependent density functional theory: 
Properties of exact functionals from an exactly-solvable model system. 



David G. Tempel' and Alan Aspuru-Guzik^ 

^Department of Physics, Harvard University, 17 Oxford Street, 02138, Cambridge, MA 
^Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford Street, 02138, Cambridge, MA 

The dissipative dynamics of many-electron systems interacting with a thermal environment has remained a 
long-standing challenge within time-dependent density functional theory (TDDFT). Recently, the formal foun- 
dations of open quantum systems time-dependent density functional theory (OQS-TDDFT) within the master 
equation approach were established. It was proven that the exact time-dependent density of a many-electron 
open quantum system evolving under a master equation can be reproduced with a closed (unitarily evolving) 
and non-interacting Kohn-Sham system. This potentially offers a great advantage over previous approaches to 
OQS-TDDFT, since with suitable functionals one could obtain the dissipative open-systems dynamics by simply 
propagating a set of Kohn-Sham orbitals as in usual TDDFT. However, the properties and exact conditions of 
such open-systems functionals are largely unknown. In the present article, we examine a simple and exactly- 
solvable model open quantum system: one electron in a harmonic well evolving under the Lindblad master 
equation. We examine two different representitive limits of the Lindblad equation (relaxation and pure dephas- 
ing) and are able to deduce a number of properties of the exact OQS-TDDFT functional. Challenges associated 
with developing approximate functionals for many-electron open quantum systems are also discussed. 



PACS numbers: 

I. INTRODUCTION 

Due to its attractive balance between accuracy and effi- 
ciency, time-dependent density functional theory (TDDFT) 
is an ideal method for computing the real-time dynamics of 
many-electron systems |T-T6|. In its original formulation by 
Runge and Gross 1 17 1, TDDFT addresses the isolated dynam- 
ics of electronic systems evolving unitarily under the time- 
dependent Schrodinger equation. However, there exist many 
situations in which the electronic degrees of freedom are not 
isolated, but must be treated as a subsystem imbedded in a 
much larger thermal bath. Several important examples in- 
clude vibrational relaxation of molecules in liquids and solid 
matrices 1 18 -20|, photo-absorption of chromophores in a pro- 
tein bath Ii21h 24I. nonlinear spectroscopy in the condensed 
phase ll25l l26l . electron-phonon coupling in single -molecule 
transport ||22l I30ti33l and exciton and energy transfer ll35l - 
[371 ED. 

In the above examples, the theory of open quantum systems 
(OQS) within the master equation approach is often used to 
model the dissipative electron dynamics ll40H43l l45l [46 1 . In 
the master equation approach, one traces over the bath de- 
grees of freedom, arriving at a simpler description in terms 
of the reduced density matrix of the electrons only. The price 
paid is that the resulting dynamics are non-unitary, and in gen- 
eral the interaction of the system with it's environment must 
be modeled in an approximate way. Even with simple system- 
bath models, the exact solution of the master equation for the 
reduced dynamics of an interacting many-electron system is 
computationally intractable. Therefore, open quantum sys- 
tems TDDFT (OQS-TDDFT) offers an attractive approach to 
the many-body open-systems problem. 

Several different formulations of OQS-TDDFT have been 
proposed in the last few years |[27ti29l [47ti5n . In Il27i. an 



OQS Runge-Gross theorem was established for Markovian 
master equations of the Lindblad form. A scheme in which the 
many-body master equation is mapped onto a non-interacting 
Kohn-Sham master equation was proposed for application to 
single-molecule transport. A Kohn-Sham master equation 
was also used in formulating the linear-response version of 
OQS-TDDFT, giving access to linewidths of environmentally 
broadened spectra f28l. In [28], OQS Casida-type equations 
were derived and used to calculate the spectrum of an atom in- 
teracting with a photon bath. A different formulation of OQS- 
TDDFT based on the stochastic Schrodinger equation rather 
than the master equation has also been developed [47 -.-.49] . 

Recently, in |[50ll5TI . the OQS Runge-Gross theorem was 
extended to arbitrary non-Markovian master equations and a 
Van Leeuwen construction was established, thereby proving 
the existence of an OQS-TDDFT Kohn-Sham scheme |52|. 
We showed that the time-dependent density of an interact- 
ing OQS can be reproduced with a non-interacting and closed 
(unitarily evolving) Kohn-Sham system. In principle, the 
closed Kohn-Sham scheme is remarkably useful for real-time 
dynamics, since it allows one to calculate any property of a 
many-body OQS by unitarily propagating a set of one-particle 
orbital equations evolving in a local potential. With suitable 
functionals, such a scheme could readily be implemented in 
existing real-time TDDFT codes HHSI. In practice, the OQS- 
TDDFT exchange-correlation potential is very complicated 
object. Not only does it have initial-state and memory de- 
pendence as in usual TDDFT 1 53 -56 1, but it must also be a 
functional of bath quantities such as the bath spectral density. 

In the present manuscript, we study exact features of the 
OQS-TDDFT closed Kohn-Sham scheme using an exactly- 
solvable one-electron model system. By focusing on a single 
electron, we are able to isolate the part of the exact functional 
arising solely from interaction with the bath, without the need 



2 



to describe electron-electron interaction effects within the sys- 
tem. 

The paper is organized as follows. In Section II, we review 
the theory of OQS for a many-electron system and the closed 
Kohn-Sham scheme presented in previous work fSOllSl |. Sec- 
tion III presents the model system to be analyzed and dis- 
cusses the procedure used to obtain exact Kohn-Sham quan- 
tities. Section IV presents results and an analysis. The paper 
concludes with an outlook and discussion of challenges for 
OQS-TDDFT in Section V. Atomic units in which e = ft = 
nie — I aie used throughout. 

n. OPEN QUANTUM SYSTEMS TDDFT USING UNITARY 
PROPAGATION 

In this section, we begin by briefly reviewing the mas- 
ter equation approach for many-body OQS. We first discuss 
the most general master equation including non-Markovian 
effects and initial correlations [62 1. We then introduce the 
Markov approximation and the widely used Lindblad master 
equation |'42'-'45', "SF, "661 . Lastly, we review the construction 
discussed in |50J, where the time-evolving density of an in- 
teracting OQS is reproduced with a noninteracting and closed 
Kohn-Sham system. 

A. The formally exact many-body master equation 

The starting point of many-body OQS is the unitary evo- 
lution for the full density matrix of the system and the reser- 
voir (we use the terms "reservoir" and "bath" interchangeably 
throughout) , 

j^p{t)^~i[H{t),pit)]- (1) 
The fuU HamUtonian is given by 

H{t)^Hs{t)+HR + V. (2) 

Here, 

Ut) = E + r-^ +Lvv.„(r„r), (3) 

1=1 i< j I ' j\ i 

is the Hamiltonian of the electronic system of interest in an 
external potential \'^x,{r,t). This potential generally consists 
of a static external potential due to the nuclei and an external 
time-dependent driving field. The system-bath coupling, V, 
is generally assumed to be weak and is treated using pertur- 
bation theory in most applications. V acts in the combined 
Hilbert space of the system and reservoir and so it couples 
the two subsystems. Hr is the Hamiltonian of the reservoir, 
which typically has a dense spectrum of eigenstates relative 



to the system. The density of states of Hr determines the de- 
cay rate of reservoir correlation functions, whose time-scale 
in turn determines the reduced system dynamics. 

Defining the reduced density operator for the electronic sys- 
tem alone by tracing over the reservoir degrees of freedom, 

Ps{t) = TrR{p{t)}, (4) 
one arrives at the formally-exact quantum master equation, 

j^psit) - -i[Hs{t),Ps{t)]+ j' dxZit - t)ps(t) +^(f). 

(5) 

Here, E(f — t) is the memory kernel and 'P(f ) arises from 
initial correlations between the system and its environment. 
It is referred to as the inhomogeneous term. The above equa- 
tion is still formally exact, as ps(f ) gives the exact expectation 
value of any observable depending only on the electronic de- 
grees of freedom. In practice, however, approximations to E 
and are required. Of particular importance in TDDFT is the 
time-dependent electronic density, 

n{r,t) = Trs{ps{t)h{r)}, (6) 

where n(r) — Y!i 8{r — f,) is the number density operator 
for the electronic system. For OQS, the continuity equation is 
not strictly satisfied and is modified to 

^«(r,f)--V.rrs{ps(f)j(r)} 

+ rr{«(r) (^^'^/TS(f-T)ps(T)+^(f))}. (7) 

Here, the first term is the divergence of the usual current 
arising in closed systems, and is referred to as the "Hamilto- 
nian current". The second term is a contribution to the current 
due to scattering of electrons with particles in the bath and 
arises from the non-unitary part of the evolution. One may 
define a "Dissipative current" by OTllTOl . 

- V • jrf„p(r,f) - Tr{n{r) {^j'^dxZ{t - x)ps{t) +m{t)^ }. 

(8) 

It will be seen later that the functional dependence of the 
Kohn-Sham potential on idispij^^t) plays an important role in 
OQS-TDDFT 

B. The Markov approximation and Lindblad master equation 

Without suitable approximations to the memory kernel 
S(f), solving eq. [sjis not easier than solving the full system- 
bath dynamics described in eq. [T] One often invokes the 
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Markov approximation, in which the memory kernel is local 
in time, i.e. 



(9) 



The Markov approximation is valid when Tj 3> is sat- 
isfied, where is the time-scale for the system to relax to 
thermal equilibrium and Tg is the longest correlation time of 
the bath. Roughly speaking, the memory of the bath can be 
neglected when describing the reduced system dynamics, be- 
cause the bath decorrelates from itself before the system has 
had a chance to evolve appreciably. Ts is inversely related to 
the magnitude of system-bath coupling, and so a weak inter- 
action between the electrons and the environment is implicit 
in this condition as well. 

The Lindblad form of the Markovian master equation. 



-Ll^L^„Ps{t)--ps{t)L, 



(10) 



is constructed to guarantee complete positivity of the den- 
sity matrix 1421 . This is desirable, since the populations of 
any physically sensible density matrix should remain positive 
during the evolution. 



As written in eq. 10 the Lindblad equation is simply a 
mathematical construction which is guaranteed to give a pos- 
itive density matrix. However, the bath operators L„„, can be 
derived microscopically starting from a system-bath Hamilto- 
nian of the form given in eq. |2]|71 1. In general, the operators 
Lnm will describe "jumps" between eigenstates of Hg induced 
by scattering of electrons with bath particles. We will discuss 
specific forms of the L„,„ in subsequent sections. 



C. Open-interacting to closed-noninteracting mapping 



Following a construction similar in spirit to f57\, one is able 
to prove that there exists a unique local potential v'(a,r,r) for 
the auxiliary system, at any interaction strength a, and for 
arbitrary E'(f,f') and ^'(f) (with some restrictions discussed 
in 1 50 1), such that 



n(r,r) = Trs{ps{t)m} = Trs{p's{t)h{r)} (13) 

is satisfied for all times. This means that for a fixed open 
and interacting system, one can always construct an auxiliary 
system with different electron-electron and system-bath in- 
teractions, such that the potential v'(a,r,f) enforces the cor- 
rect density evolution. If one sets a = 0, but keeps E'(f) — 
S(f) and ^'(f) = ^(f) then the auxiliary system is a non- 
interacting, but open Kohn-Sham system. This is similar to 
the construction used in ll27l l28l . but encompasses the non- 
Markovian case as well. However, one may also choose a = 
and S'(f) = ^'{t) — 0, whereby the density of the original 
open system is reproduced with a closed and non-interacting 
Kohn-Sham system. In this case, the density matrix in the 
auxiliary system is a pure state given by 



p^(f)=4>*(f)4>(0, 



(14) 



where <I>(f ) is a single Slater determinant. This determinant 
is constructed by propagating a set of single-particle orbital 
equations. 



i^0,(r,f)|-^v2 + v,,(r,f)|0,(r,O (15) 

as in usual TDDFT. The density of the original interacting 
and open system is then simply obtained by square-summing 
the orbitals. i.e. 



In this section, we briefly outline the proof given in |30l, 
whereby the density of an interacting and open electronic sys- 
tem is reproduced using a closed and non-interacting Kohn- 
Sham system. 

Starting from the many-body master equation given in eq.[5] 
one considers an auxiliary "primed" system described by a 
density matrix p'g{t), evolving with a different Hamiltonian 

^K0--^Ev2 + f ^^+^v'(a,r„0, (11) 



i=l 



at interaction strength a, and with different memory ker- 
nel S'(f) and initial correlations ^'(t). This auxiliary system 
evolves under the master equation 



|p^(f) = -.[7f^(f),p^(f)]+£«'TS'(f-T)p^(T)+>l"(f). 



n{r,t) ^Y.\'^i{r,t)\^ ^Trs{ps{t)h(r)). (16) 

In analogy to usual TDDFT, the Kohn-Sham potential is 
partitioned as 



vfa (r, f ) = v„, (r , f ) + V,, (r , f ) + (r, f ) , 



(17) 



(12) 



where v/,(r,f) is the Hartree potential and the unknown 
functional v"c^"{r,t) accounts for electron-electron interac- 
tion within the system as well as interaction between the sys- 
tem and bath. Ideally, these two contributions well be addi- 
tive, so that one may use standard adiabatic functionals to ac- 
count for electron-electron interaction within the system and 
construct dissipative bath functionals to account for system- 
bath coupling. In general. 



v'ir{r,t) = v;r"(r,f)[«,S,^,ps(0),4>(0)]. (18) 
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Formally, the open-systems exchange-correlation potential 
is a functional not only of the density, but also of the memory 
kernel, inhomogeneous term and initial state of the interacting 
open system, as well as the initial state of the closed Kohn- 
Sham system. 



m. AN EXACTLY-SOLVABLE MODEL SYSTEM 

In this section, we construct the exact OQS-TDDFT Kohn- 
Sham potential for an exactly-solvable model system: one 
electron in a harmonic well evolving under the Lindblad equa- 
tion. Our analysis focuses on two limiting cases of the Lind- 
blad master equation. The first limit is that of pure dephasing 
without relaxation in which the bath decoheres the system, but 
no energy is exchanged. The second limit is that of relaxation 
with no pure dephasing. 

A. Construction of the exact OQS functional 

With the system Hamiltoiuan given by 



the Lindblad equation 



1 

2 dx^ 



1 



2 2 
(OX, 



(19) 



-i[Hs,Ps{t)] 

fim I 



(20) 



n{x,t) = \^{x,t)\^ ^Trs{ps{t)^{^)}- 



(23) 



Here, 



vtv(x,0 = veAx)+n{x,t) + vT"{x,t) (24) 

is the OQS-TDDFT Kohn-Sham potential. We now sub- 



stitute eq. 21 in eq. 22 and use the fact that for one electron. 



the exchange potential exactly cancels the self interaction in 
the Hartree potential. The result is an exact expression for the 
OQS-TDDFT correlation potential in terms of known quanti- 
ties. 



open 



{x,t) 



d_ 
dt 
1 



a{x,t) 



dx 



a{x,t) 



4n{x,t) dx^ 



1 



in{x,ty 



n{x,t) 



-n{x,t) 



Vextix). (25) 



As in Il57ll59l , it is instructive to separate v"''™(x,f ) into an 
adiabatic and dynamical part. 



y«P™(x,f)=vf"(x,f)+vf (X,f) 



Here, 



vf{x,t)- 



An{x,t) dx^ 



n(x,t) - 



8n{x,ty 



-n{x,t) 



(26) 



(27) 



is the exact functional of ground-state DFT for one elec- 
tron, or two electrons in a spin-singlet evaluated on the in- 
stantaneous OQS density 1611 . The dynamical part. 



can be solved exactly to obtain ps{t)- The exact time- 
dependent OQS density can then be constructed using eq. [6] 
With an exact OQS density, it is a simple exercise to construct 
the closed Kohn-Sham system which reproduces this density 
using a unitary evolution. In this section we work in one di- 
mension, but the formulas apply to higher dimensions as well. 

For one electron, there is a single occupied Kohn-Sham or- 
bital given by 



(21) 



By construction, this orbital must evolve under the time- 
dependent Kohn-Sham equation 



dt 



1 d^ 

4>{x,t) - { ---j-2+v,{x,t)}^{x,t), (22) 



in such a way that the true OQS density is reproduced for 
all times, i.e. 



viy"ixj) = -^a{x,t) 



-a{x,t) 



(28) 



is a strictly-dynamical contribution which vanishes when 
the system is in thermal equilibrium. From the generalized 
continuity equation, eq. [7] one finds that the phase of the 
Kohn-Sham orbital is given by 



dx n[x,t) n[x,t) 



(29) 



We see that for OQS, the dynamical correlation potential 
vp''{x,t) contains a contribution from the Hamiltonian cur- 
rent that is also present in usual TDDFT. In addition, there is 
a new contribution from the bath arising through j^,-^.^^. There- 
fore, the OQS correlation potential is a functional not only of 
the Hamiltonian current, but the dissipative current as well. In 
the following, we will consider the two limiting cases men- 
tioned above: Pure dephasing without relaxation and relax- 
ation without pure dephasing. 
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B. Pure dephasing without relaxation 



The coherences evolve according to 



We first consider a situation in which the OQS evolves un- 
der a Lindblad master equation which induces pure dephasing, 
but no relaxation 138] |^ |63l |64l . In this case, the Lindblad 



operators are diagonal in a basis of eigenstates of eq. 19 and 
given by 



Ami = Snm\l y k^X'"!, 



(30) 



where \m) is the mth eigenstate of the oscillator. Since the 
operators L„, are diagonal, the populations Pmm{t) remain un- 
changed as the system evolves. This implies that energy is 
conserved during the evolution, as the bath cannot drain en- 
ergy away from the system. However, the coherences given 
by the off-diagonal density matrix elements, p„m{t), decay ex- 
ponentially with a rate given by 



'^decoherence ~ 2^"^'"^"^"^' 



Pure dephasing describes a situation in which system-bath 
collisions are elastic, so that the bath decoheres the system 
without exchanging energy. 



C. Relaxation without pure dephasing 

Next, we consider the case of relaxation without pure de- 
phasing, in which the Lindblad operators are strictly off- 
diagonal ll64l . i.e. 



L,nn = \f%^i\m) (« I , m^n 



(32) 



In order to ensure that the populations of the equilibrium 
solution obey detailed balance, one also requires that 



^™=e'^'"'""7™, (33) 

where j3 = is the inverse temperature. With the Lind- 
blad operators given by eq. |32] the populations evolve accord- 
ing to 



dt 



Pnn{t) ^Y^Ynme Pmm{t) - P„n{t)Y,7mne . (34) 



The first term in eq. 34 expresses the rate at which popu- 
lation is transferred to p„„ from all other populations m ^ n, 
while the second term expresses the rate at which population 
leaves Pjj,,. It can be readily verified that the right hand side 



of eq. 



34 



vanishes at equilibrium where p, 



eq 



these two rates balance. 



and 



dt 



Pnm{t) = { - 



/ 



it)- 

(35) 

This shows that even in the absence of pure dephasing, re- 
laxation still necessarily implies decoherence. For the special 
case of a two-level system, one can readily see from eq's. 34 



and 35 that the decoherence rate is exactly half of the relax- 
ation rate. This is analogous to the ubiquitous phenomenolog- 
ical formula from NMR spectroscopy 



1 

7^ 



1 



^2 



(36) 



relating the time-scale for decoherence T2, to that of relax- 
ation Ti and pure dephasing T2 ll64l l65l . In the absence of 
pure dephasing, the time-scale for decoherence is twice that 



(31) of relaxation, which is confirmed by eq's. 34 and 35 



IV. RESULTS AND ANALYSIS 

We now present and analyze the results of the inversion pro- 
cedure mentioned in the previous section for obtaining the ex- 
act Kohn-Sham quantities. For all calculations, we choose the 
initial state of the OQS to be the pure state 



\wm 



:(|0> + |1)), 



(37) 



which corresponds to a density matrix with initial elements 
Poo(O) = Pii(O) = Pio(O) = Poi(O) = 5 and all other entries 
equal to zero. The frequency of the oscillator is taken to be 
co = l. 

For the pure dephasing case (section III B.), we choose the 
parameters 7o = 7i = 0. 15 a.u., leading to a decay of the initial 



decoherence 



0.15 



coherence between |0) and |1) at a rate of T^' 
a.u. as given by eq.[3T| 

For the case of relaxation without pure dephasing (section 
III C), we choose the population transfer rate from |1) to 
|0) to be 7oi =0.3 a.u. and the inverse temperature to be 
ji — I- Using eq. 33 this gives 7io = ^ ~ 0.11 a.u. and 
720 = ^ ~ 0.04 a.u. At this relatively low temperature, the 
transfer rates to excited states higher than |2) are sufficiently 
small that they can be neglected. Using eq 
of the initial coherence between |0) and |1) 



35 the decay 



decoherence 



_ roi 



is found to be 
2.1 a.u. These parameters all cor- 



respond to underdamped motion in which the decay to equi- 
librium occurs on a much longer time-scale than the oscilla- 
tion period of the system. This condition is implicit in the 
assumption of weak system-bath coupling. 
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A. The density, current and dissipative current 

In Figures [T] and |2] we plot the density, current density 
and dissipative current density as a function of time for the 
pure dephasing and relaxation without pure dephasing mas- 
ter equations respectively. Since we have chosen parameters 
corresponding to weak system-bath coupling, the dissipative 
current j^isp {x,t)is significantly smaller than the Hamiltonian 
cun-ent j{x,t), by approximately a factor of if^atherence 
the pure dephasing case. In the figures, we have multiplied 
jdisp{x,t) by a factor of 5 to make it more visible. 

In Figure[T] we see that the coherences decay exponentially 
in time as the density oscillates in the harmonic well, while the 
populations remain unchanged, as expected for pure dephas- 
ing. In the absence of coupling to the bath, time translational 
invariance would imply that the snapshots in the right hand 
column would be the same as those on the left, which occur at 
a time t = 2n earlier. The loss of coherence manifests itself in 
the decay of the currents and changing density profile as the 
system evolves toward the fully mixed state poo = pn — ^, 
Pio — Poi = at equilibrium, which is an incoherent sum of 
the ground and first excited states. One can also see that the 
dissipative current is proportional to the spatial integral of the 
real part of the coherence, and lags the hamiltonian current by 
a phase of | which is proportional the spatial integral of the 
imaginary part. Only the real part of the coherence contributes 
to the density evolution. 

The contribution to the density from the coherences is an- 
tisymmetric about the origin, while the contribution from the 
populations is symmetric. When these two contributions are 
superposed, the density acquires an asymmetric profile, which 
can be seen at integer multiples of n where the coherent con- 
tribution is a maximum. At half integer multiples of n, the 
contribution from the real part of the coherences instanta- 
neously vanishes and the density becomes symmetric. In con- 
trast, the currents are always perfectly symmetric about the 
origin since they are proportional only to the spatial integral 
of the antisymmetric coherences. 

In Figure |2] we see that in addition to decoherence there is 
population transfer, since we have included energy relaxation 
in the master equation. By t — the system has already 
nearly approached the equilibrium state 

1 e-l^'^ 
PsC = -) = 10) (01 + Y^e^ 1 1) (1 1 ' (38) 

which for the parameters we have chosen corresponds to 
Poo(f = oo) « 0.73 and P\\{t = °°) ~ 0.27. Since the equi- 
librium density is dominated by the ground-state, it is nearly 
gaussian, but slightly flattened due to a mixing in of the first 
excited state. The dissipative current has a similar profile to 
that seen in li31J . which results from the fact that population 
transfer generates no Hamitonian current, but does give rise 
to a change in the density which must be compensated by 

jdi.sp (-^lO- 



Figure 1: Time evolution of the density n{x,t) (red-solid), the Hamil- 
tonian current j{x,t) (green-dot-dashed) and the dissipative current 
jdisp(x,t) (blue-dotted), for the pure dephasing master equation. 
jdisp(x,t) has been scaled by a factor of 5 to make it more visible 
in the figure. Vgxt (x) = j (O^x^ (green-dashed) is shown for reference 
as well. 



B. The exact Kohn-Sham potential, Kohn-Sham orbitals and 
density matrices 

In Figures |3] and |4j we show the full Kohn-Sham potential 
and its dynamical contribution (eq. [28] l for the density evolu- 
tion presented in Figures[T]and|2]respectively. In Figure|5] we 
plot the density and dynamical potential vp'"{x,t) as a func- 
tion of space and time for the pure dephasing case as well. 
Examining Figures [3] and |5] we see that for pure dephasing, 
the dynamical potential becomes extremely repulsive in re- 
gions where the contribution to the density arising from the 
coherences in the OQS density matrix is large. For instance, 
at f = ;r, the contribution to the density from the coherences is 
positive for X < and negative for x > 0. Likewise, vp"{x,t) 
is positive forx < and negative forx > 0, which corresponds 
to a field which will tend to suppress the coherent part of the 
density. For t — 2n, the sign of the coherent part of the den- 
sity is reversed and vp"{x,t) is as well, again corresponding 
to a potential which redistributes the density so as to suppress 
its coherent contribution. This behavior gives us some insight 
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Figure 2: Time evolution of the density n{x,t) (red-solid), the Hamil- 
tonian current j{x,t) (green-dot-dashed) and the dissipative current 
jdisp{^,t) (blue-dotted), for the master equation including relaxation 
with no pure dephasing. jdispi^^t) has been scaled by a factor of 5 to 
make it more visible in the figure, v^.^, green-dashed) 
is shown for reference as well. 



into how the Kohn-Sham system is able to use a unitary evolu- 
tion to mimic the effect of the bath. In the Kohn-Sham system, 
decoherence is effectively converted into a control problem, in 
which the local time-dependent field vi^"{x,t) drives the den- 
sity in the same way that collisions with the bath would in 
the true open system. This is a complicated task, since not 
only must Vc^"{x, t) drive down the coherences, but it must do 
so without affecting the contribution to the density from the 
populations. The situation is somewhat different in Figure |4] 
since now the potential not only suppresses the coherent den- 
sity evolution, but it also transfers population from the first 
excited state to the groundstate to reach the thermal equilib- 
rium density distribution. 

At f = 0, the Kohn-Sham potentials appear to be extremely 
large, however this occurs in a region of space where the den- 
sity vanishes and has little effect on the dynamics. At long 
times, one sees that the dynamical potential decays, while in 
the pure dephasing case, the full Kohn-Sham potential ac- 
quires a double well structure. This corresponds to the adi- 



mixed, equilibrium-state density. In figure]?] the double well 
structure is less pronounced, since the equilibrium state is 
dominated by the groundstate for which the adiabatic poten- 
tial reduces to Vg^t (up to a constant), which is parabolic. 
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Figure 3: Time evolution of the density n(x,t) (red-solid), the 
full Kohn-Sham potential Vj(,s(x,?) = Vc'''™(x,?) + Vexi{x) (green-dot- 
dashed) and the dynamical correlation potential vi^"{x,t) (blue- 
dotted) for the pure dephasing master equation. Note the change of 
scale for the t = frame, v^y (x) 
for reference as well. 



^ co^x^ (green-dashed) is shown 



abatic Kohn-Sham potential (eq. 27 1 evaluated on the fully 



In Figures |6] and |7j the Kohn-Sham orbital corresponding 
to the two different master equations is shown. The imagi- 
nary part of the orbital arises from the Hamitonian current as 
in usual TDDFT, but also has a contribution from the dissi- 
pative current. At long times, the currents decay and so the 
imaginary part vanishes while the real part is simply given by 
^/njxj). In the long-time limit, this orbital must reduce to 
an eigenstate of the adiabatic Kohn-Sham potential evaluated 
on the equilibrium OQS density. In this way, the Kohn-Sham 
system reproduces the coiTect equilibrium density of a mixed- 
state density matrix using a pure-state wavefunction. In prin- 
ciple, this wavefunction need not be the groundstate, but it is 
in the cases studied here. 

Lastly, we show the real pait of the real-space density ma- 
trices for both the Kohn-Sham system and the true OQS in 
Figures |8] and |9] In the Kohn-Sham case, this is simply the 
pure-state density matrix 
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Figure 4: Time evolution of tlie density n{x,t) (red-solid), the 
full Kohn-Sham potential v^.j(x, ?) = v°'""^{x,t) + Vext{^) (green-dot- 
dashed) and the dynamical correlation potential vi^'"{x,t) (blue- 
dotted) for the master equation including relaxation with no pure 
dephasing. Note the change of scale for the f = frame. Vext{x) = 
green-dashed) is shown for reference as well. 



j dx j dx'\pks{x,x ,t)\^ = 1, (41) 
while for the true system, 

j dx j dx'\popen{x,x ,t)\^ < 1. (42) 

The Kohn-Sham density matrix is able to satisfy eq. |4T|by 
increasing the off-diagonal regions and thus compensating for 
the lost volume in the true density matrix without affecting the 
density. 

C. The Markovian bath functional 

Our analysis of the exact OQS-TDDFT functional in sec- 
tion IV B also sheds some light on the Markovian bath func- 
tional (MBF) that we first presented in |50|. In the MBF, one 
approximates the dynamical correlation potential as 

" {x,t)^ vf "'^^^ (x, t)^X ^j^^dx', (43) 

where A is a fitted parameter. In this approximation, one 
neglects the contiibution to vv "(x,f) arising from the second 
term in eq. |28] involving gradients of the phase and one also 
neglects the contribution from the dissipative current. The first 
approximation is valid when the phases are spatially slowly 
varying, while the second approximation is valid if the cou- 
pling to the bath is weak. Lastly, one assumes 



PkAxy,t)^rix,t)<pix't), (39) 

while in the OQS case it is the position representation of 
the solution to eq. [20l 

Popen{x,x' ,t) = {x\ps{t)\x') . (40) 

Since we start in a pure-state at f = 0, the Kohn-Sham den- 
sity matrix coincides with that of the true OQS. As the sys- 
tem evolves, the two begin to differ as the true density ma- 
trix looses purity while that of the Kohn-Sham system does 
not. However, it can be seen that at all times piis{x,x,t) = 
Popen{x,x,t) = n{x,t) and the Kohn-Sham system reproduces 
the true OQS density on its diagonal. 

In both Figure [8] and Figure |9] we see that the Kohn-Sham 
density matrix has large off-diagonal support which is not 
present in the true density matrix. This is consistent with the 
fact that the Kohn-Sham density matrix must remain idempo- 
tent. i.e. 



—a{x,t)o^-Xa{x,t). (44) 

i.e. the phase of the orbital is exponentially damped, which 
is valid in the Markovian limit. The MBF is a simple yet prac- 
tical functional as it only involves knowledge of the Hamilto- 
nian cuiTent and density at each instant in time. It can readily 
be implemented in real-time codes with existing propagation 
schemes [78|. Numerical simulations of dissipation of ex- 
cited states of a model Helium atom show promising results, 
notwithstanding the simplicity of the MBF ll50l . 

In Figure |5] we show a comparison of the MBF 
vp"'^^^ {x,t) and the exact functional v'P"{x,t) for the pure 
dephasing master equation of the model system presented in 
the previous sections. The fitted decay parameter is taken 
to be the decoherence rate in the master equation, i.e. A — 
'^decoherence =0-15. This is natural since it is the only decoher- 
ence timescale in the problem. The MBF is seen to reproduce 
the long-time behavior of the exact dissipative potential rea- 
sonably well, however it deviates significantly for short times. 
This arises from the fact that the Markov approximation only 
provides an adequate description of the dynamics for times 
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Figure 6: Time evolution of the density n{x,t) (red-solid), the real 
part of the Kohn-Sham orbital 1ie(j){x,t) (green-dot-dashed) and the 
imaginary part of the Kohn-Sham orbital 3m<j){x,t) (blue-dotted), for 
the pure dephasing master equation. Vext (x) = j oP'x^ (green-dashed) 
is shown for reference as well. 
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Figure 7: Time evolution of the density n{x,t) (red-solid), the real 
part of the Kohn-Sham orbital 1ie(j){x,t) (green-dot-dashed) and the 
imaginary part of the Kohn-Sham orbital 3m(j){x,t) (blue-dotted), 
for the master equation including relaxation with no pure dephasing. 
Vext (x) = ifflV ( green-dashed) is shown for reference as well. 
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Figure 8: Time evolution of the real part of the Kohn-Sham den- 
sity matrix 9lep^.j(x,x',f) (left column) and true OQS density matrix 
5Repop(,„(x,x',f) (right column) for the pure dephasing master equa- 
tion. 



Kohn— Sham True 




Figure 9: Time evolution of the real part of the Kohn-Sham den- 
sity matrix 9lep^.j(x,y ,f) (left column) and true OQS density matrix 
'3iepopg„{x,x' ,t) (right column) for the master equation including re- 
laxation with no pure dephasing. 



longer than the decoherence time. At much shorter times it is 
not in general valid. In Figure [5] we see that this is roughly 
the timescale for the MBF to become accurate. Also, because 
the MBF neglects the dissipative current which is out of phase 
with the Hamiltonian current, there is a slight phase shift in 
the MBF with respect to the exact dissipative functional which 
contains both current contributions. 



V. CONCLUSION AND OUTLOOK 

In summary, we have explored the behaviour of the 
exchange-correlation potential for dissipative open quantum 
systems using an exactly-solvable one-electron system. The 
two limiting cases (pure dephasing and relaxation) provided 
insights into the time-dependence and form of the dissipa- 
tive potential that will need to be described in OQS-TDDFT 
functionals. We explored the behavior of the closed auxil- 
iary Kohn-Sham system, as the system evolves under decoher- 
ence and relaxation. This is valuable information for the de- 
velopment of many-body, realistic bath functionals for OQS- 
TDDFT. 

We have seen that the Kohn-Sham system mimics the ef- 



fect of a bath by unitarily evolving with a dynamical poten- 
tial, which depends on both the Hamitonian and dissipative 
cuiTents. For a pure dephasing master equation, this potential 
tends to suppress the coherent part of the density evolution in 
a delicate way so that energy is still conserved and popula- 
tions remain unchanged. When relaxation is also present, the 
potential drains energy away from the system and reaches an 
equilibrium distribution which satisfies detailed balance. 

When dealing with a many-electron system, one expects 
that the functionals should retain some of the same features 
as in the one-electron case studied here, however, there are 
some fundamental differences. Firstly, one does not know 
the many-electron eigenstates a priori and so it is generally 
not possible to explicitly write down Lindblad operators in 
terms of these eigenstates as was done in eq. 20 Instead, 



one should start from a master equation written in terms of 
effective single-particle eigenstates, such as a basis of or- 
bitals from an equilibrium-state Kohn-Sham-Mermin calcu- 
lation [27,18., 72 1 . 

Secondly, the simple form of the dynamical functional ex- 



pressed in eq. 28 no longer holds. In the simulation we have 
considered here, there are only two relevant levels involved in 
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the dynamics. In general, for a many-level system the bath 
will induce different relaxation and dephasing rates for dif- 
ferent eigenfrequencies of the system, depending on the bath 
spectral density Il42l [44l [64 \ . Rather than simply coupling to 



the currents in a semi-time-local manner as in eq. 28 the ex 



act functional will need to damp different Fourier components 
of the current at different rates. In the time-domian, this cor- 
responds to a complicated memory-dependence on both the 
Hamiltonian and dissipative currents at earlier times. The 
form of this memory-dependence should clearly depend on 
the bath spectral density, but its exact structure will need to be 
investigated in future work. 

In some respects, OQS-TDDFT functionals are similar to 
existing current-dependent functionals in TDCDFT, where 
frequency-dependent dissipation arises due to coupling to cur- 
rents at earlier times via a stress tensor ||73H771 . However, 
the physical origin of the dissipation is very different in the 
two theories. In TDCDFT, the frequency dependence of the 
stress tensor depends on viscoelastic coefficients derived from 
the uniform electron gas. In contrast, the OQS-TDDFT func- 
tional depends on the spectral density of a bosonic bath such 
as phonons or photons. 
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